Random loose packing and an order parameter for the parking lot model 
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Abstract 

We have obtained the random loose packing fraction of the parking lot model (PLM) by taking 

the limit of infinite compactivity in the two- variable statistical description of Tarjus and Viot for the 

PLM. The PLM is a stochastic model of adsorption and desorption of particles on a substrate that 

have been used as a model for compaction of granular materials. An order parameter p is introduced 

^^ , to characterize how far from a steady state situation the model is. Thus, configurations with p < 1 

' age. We propose that p can be a starting point in order to stablish a connection between Edwards' 

^~>| statistical mechanics and granular hydrodynamics. 
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1 Introduction 



More than twenty years ago Edwards and Oakeshott proposed a statistical mechanics framework for 

granular materials in mechanical equilibrium jlj . The idea was to replace the energy by the total volume 

f"j^ ' of the sample V. Thus the entropy of the system is defined as 5 = log fi, where ft is the number of stable 

O . states for a given volume V (and given number of grains N). The quantity equivalent to temperature 

within this description is the compactivity x = {dS/dV)~^ . In connection with Edwards proposal, Aste 

j^ . et al has found for static packings of spheres an invariant distribution of Voronoi volumes at the grain 

level [5]. Also, the random close packing fraction (f>rcp and the random loose packing fraction (prip have 

been associated with configurations with x ~^ and x — ^ oo, respectively [3j|4j[5]. 

' O ' It has been shown recently that two granular samples with the same packing fraction may not have the 

^ , same properties, and an aditional macroscopic variable must be introduced in the statistical description. 

In reference [B] was found that we can produce two samples with the same but different stresses. Indeed, 

the stress tensor a has been included in a more general statistical description [7] (see [B] and references 

therein). 

We deal in this article with the so called parking lot model (PLM) introduced by Nowak et al in 
iy-\ ' the context of experiments of compaction of grains [8]. Starting with the statistical mechanics for the 

PLM proposed by Tarjus and Viot j9], we obtain the random loose packing fraction for this model and 
^»0 ' introduce an order parameter p that characterizes how far from a steady state situation the model is. We 

propose that a quantity analog to p can be used as the order parameter in the continuum description of 
lO ' slow and dense granular flows by Aranson and Tsimring, mediating how solid and fluid is the form of the 

^^ , stress tensor [TUl [TT] . 

The PLM is a model of random adsorption and desorption of particles on a substrate. Particles are 
disorbed with rate p- and adsorbed with rate p^, with a no overlapping condition. For a given initial 
condition of the substrate, the model converges to a stationary state packing fraction 0e around which it 
fluctuates. 0e depends only on the ratio K = p^/p^, which allow us to map the parameter K ioT or Q 
\^ • of references [3 [12]. For large K we have (j)e = 1 — 1/ log iiT [13] and the convergence to a stationary state 

5^ , is very slow, reminiscent of glassy behavior [T^ . The PLM is though to represent an average column of 

grains. 

Tarjus and Viot characterized the configurations produced by the PLM with two variables I9j. One 
of them is the packing fraction (j> and the other is the insertion probability $, which is the available 
line fraction for a new insertion. Tarjus and Viot recognized the need for an additional variable because 
of some memory effects observed in experiments that can be reproduced within the PLM if we change 
K in the course of a Monte Carlo simulation: we can generate two conflgurations with the same but 
with different subsequent evolution and by allowing a variation in K, for a given finite time, we can 
obtain more dense substrates [T5]. Thus, the history of the configuration is encoded in $, which is a 
structure variable |16) . The aditional variable $ is only needed to characterize configurations which are 
not produced in steady state, since in steady state the insertion probability $e is given by (eq. 2 in [16): 

$e = (l-0)exph0/(l -(/.)]. (1) 
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Figure 1: Distribution of Voronoi lenghts of the PLM for K^ip ~ 20 obtained with Monte Carlo Sim- 
ulations (circles) . The continuum line is a gamma distribution with shape parameter A: = 2 [21 [20] . In 
the inset we show 2/3 — 0r/p vs. the insertion probability $ for the PLM as given by eqs. [7|and[Hl The 
dotted line is equation [3| At $ « 0.01 the difference in i^rip is less than 1% . 

2 Results 

If A is the total lenght available for a new insertion (A = $1/, where L is the size of the system), for 
given iV, L and A we have a configurational integral in terms of the gaps hi (eq. 25 in [S]): 



Z = 



f ... /" dh^5lL-N-f2h^]slA-f2^ih,~l)ih,-l)\ 



(2) 



with 9 being the step function. Equation [5| can be solved with a saddle point method in the limit of large 
iV, L and A, with (f) = N/L and $ = A/L fixed. 

In steady state, for K ~ 100 we have $ k, 0.01 15 . In the limit of small $ we obtain a random loose 
packing fraction by taking the limit x ~^ cxd in the above description (see Appendix): 



4>rlp 



$1/2^ 



(3) 



Since the random loose packing of granular media depends on friction [17j , from equation [3| we see 
that greater <I> suggest grains with larger friction coefficient. For a given 0, a greater $ can also be 
associated to greater heterogeneity of voids at the grain level. In the limit $ ^^ 0, from equation[3|we get 
(j)rip = 'i-l'i- This value for (j)rip correspond to smooth grains. A gamma distribution of Voronoi lenghts 
[20] is obtained for K^ip ~ 20 (see figure [1]). The invariant distribution found in experiments by Aste et 
al for spheres packings in mechanical equilibrium is also a gamma distribution [2] . 

It have been reported that at a packing fraction near the 4>rip found by us the process of compaction 
is slowered. Before reaching (/)e, typically we have four differents regimes as we record the evolution of 
the packing fraction in a Monte Carlo simulation of the PLM [18]. During a first stage, (j) "... increases 
rapidly until a value of around 0.65" [18], and from this point afterwards the increase in is considerably 
slower. Thus, (prip can play a role in the onset of jamming in granular materials, as have been speculated 
in reference [l9] . 

We introduce now a parameter p which is, basically, the quotient between Z given by equation [5] and 
Zt which is the configurational integral without the restriction of having a definite value of A, i.e. when 
only the first 5 in equation [2] is considered. This leads to: 



p(0,$) =cxp(s -St), 



(4) 




0.84 



0.82 



B0.78 



0.76 



0.74 



0.7 0.75 0.8 0.85 0.9 



0.72 




0.005 0.01 0.015 0.02 



Figure 2: a) Parameter p (eq. [3]) as a function of packing fraction (j) for different values of the insertion 
probability $. From right to left: $ = 0,0.001,0.002, ...,0.018,0.019. b) 1-p on a logarithmic scale as a 
function of for different values of $ (same data than in a)). From the intersection with the horizontal 
axis we can estimate the maximum value of packing fraction (pmax for a given $. c) Circles: ipmax vs. $ 
as obtained from b). When (f) -^ (f>max, then p ^>- 1. The solid line is equation [U the steady state relation 
between (j) and $. 

with St = (p + (p\n ( — ^ ] [S] and the entropy density s is given by equation [5] in the Appendix. In figure 

[5^ we can see p{4>, $). For a given value of $, p — ?> 1 as is increased up to a maximum value 4>max- 
This maximum value of </> can be estimated from a graph like the one shown in figure [2b- 

In figure [2}: it can be seen that p — )> 1 is equivalent to say that we are approaching a steady state 
situation. With this in mind, in figure [3] we plot the insertion probability $ vs. the packing fraction </> 
for a given value of the parameter p. We did this by solving numerically the relevant equations needed 
for to evaluate equation H) It is worth to remember that a (</>, $) statistical description for the PLM 
makes sense only if we need to consider configurations out of steady state, since in steady state these two 
variables are related by equation [TJ Thus, configurations with p < 1 age. 



3 Discussion 

In the last stages of a Monte Carlo simulation of the PLM is increasing into (pe with a small variation 
in the insertion probability $. Talbot, Tarjus and Viot found in reference [TS] that for K > Kc « 100, 
there is a minimum in $ as a function of time that occurs at a packing fraction 0™ (K) . This implies that 
for slow compaction ii K > K^ the system can increase its packing fraction while increasing or decreasing 
$, depending if </> is greater or smaller than </>„• This should have consequences for processes (curves) 
in a (0, $) plane like the one shown in figure [31 Thus, for > (p^-ip in the PLM we have two zones in 
which we expect different behaviour. In terms of K these zones are: Krip < K < Kc and K > Kc, with 
Krip ~ 20 and Kc ~ 100. 

In reference |21j the authors reported a phase transition when inserting slowly a rod into a column 
of grains: the system's response to shear changes at a certain packing fraction ipc- ipc can be localized 
by monitoring the change in height Ah of the column, after removing the rod, as a function of </> [21] . It 
would be interesting to put the results of reference [H] in terms of a process that starts from a packing 
fraction (pi on the p = 1 curve of figure [3] and ends on a packing fraction (j)f, with (pi > cpf |21) . 

Can this two- variable description of the PLM be related to the volume-stress proposal of Edwards 
and others? We believe that a connection can be made by using a quantity analog to p{(f>, $) given by 
equation [4| as the order parameter in the stress tensor aij proposed by Aranson and Tsimring in their 
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Figure 3: Insertion probability $ vs packing fraction cj) for different values of the parameter p. The 
particular case p = 1 is given by equation [H the steady state relation between and $. The dark region 
are not available (0, $) configurations and its frontier is defined by the p = 1 curve. Configurations with 
p < I age. The dotted line is a schematic representation of 4>m(K) [13]. 



continuum description of slow and dense granular flows |I0| lllj : 



M 






dxi 



+ a%[p+il-p)S., 



(5) 



where c" is the stress under static conditions with the same geometry. Thus, p in eq. |4] can control 
how fluid and solid is the form for the stress tensor. For p < I we have a partially fluidized granular 
medium. At best, this is a first step, a suggestion, towards a real connection between Edwards' statistical 
mechanics and granular hydrodynamics 22 . 

Finally, since configurations with p < I in the PLM age we must say something about the relevant 
time scales for this model. Kolan, Nowak and Tkachenko have found that the low relaxation frequency ujl 
and the high relaxation frequency ujh for this model are given by ui^ — -p- exp(— 2/(5e) and ujh = p+S, 
with S — {1 — (j))/4> [11]. Thus, in order to speak of thermodynamic processes in figure 121 the observation 
time tabs [23] must satisfy 1/ujh ^ tobs ^ l/^^L- From figure 5 of reference [13], we can see that this 
condition on tobs can be satisfied only for high packing fractions. Only for K > Kc we have at least two 
orders of magnitude of separation between ujl and ujh- 

4 Conclusion 



We have obtained in this article the random loose packing fraction (j)rip for the parking lot model, where 
^rip is the lower packing fraction in which we can find a sample in mechanical equilibrium. We have done 
this by taking the limit of infinite compactivity in the statistical description of Tarjus and Viot, in which a 
macro state is characterized by its packing fraction (j) and its insertion probability $. The compactivity x 
is the analog of temperature in the statistical mechanics for granular materials proposed by Edwards. We 
have proposed an order parameter p{(f), $) that characterizes how far from a steady situation the model 
is. Thus, configurations with p < 1 age. With p, we proposed a connection of statistical mechanics with 
the continuum description of slow and dense granular fiows by Aranson and Tsimring. By considering 
the relevant time scales for this model obtained by Kolan, Nowak and Tkachenko, we have argued that 
for blocked configurations (0 > (prip), only for even higher packing fractions we expect to be able to speak 
of thermodynamic processes for this model. 
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Figure 4: Equations 29 and 30 of reference j9]. For a given insertion probability $, we have a maximum 
z, which imphes a maximum packing fraction <^max- The horizontal line corresponds to = 2/3, the 
maximum random loose packing fraction for the PLM (see eq. [3]). 
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6 Appendix 



Tarjus and Viot obtained the entropy density s{(j), $) (Z ~ exp(Ls)) (eq. 28 in [S]): 

' z + y[l — cxp(— z)] ' 



(1 - (f))z + y^ + (j)ln 



z{z + y) 



(6) 



with z — z{(j), $) and y = y{(j), $) being solutions to the coupled equations 29 and 30 of reference [5], which 
can be seen in figure ID We have that z = {dS/dL)N^A can be interpreted as the inverse of compactivity 
and y — {dS/dA)N,L- 

We can obtain the random loose packing fraction of this model by considering the limit z — > (x — >■ cx)) 
[3j m [5] in equations 29 and 30 of reference [9] . We obtain: 



$ = 





1 


(i + y) 


1 1 71 1 "^ 


[^ ' y ' 2(1+!/) J 



and 



= 1 



(l^rip ' y 2(1 + 2/) 

from which we can eliminate y to obtain (/)f./p($), as can be seen in the inset of figure [T] 
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